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Abstract 

(N 

Random flights in R d ,d > 2, with Dirichlet-distributed displacements and uniformly 

distributed orientation are analyzed. The explicit characteristic functions of the position 

X d (t), t > 0, when the number of changes of direction is fixed are obtained. The probability 

Qv distributions are derived by inverting the characteristic functions for all dimensions d of R d 

£SJ and many properties of the probabilistic structure of X d (t), t > 0, are examined. 

If the number of changes of direction is randomized by means of a fractional Poisson 
process, we are able to obtain explicit distributions for P{X d (t) £ dx^} f° r au d > 2. A 
Section is devoted to random flights in R 3 where the general results are discussed. 
Ph The existing literature is compared with the results of this paper where in our view 

the classical Pearson's problem of random flights is resolved by suitably randomizing the 
step lengths. The random flights where changes of direction are governed by a homoge- 
neous Poisson process are analyzed and compared with the model of Dirichlet-distributed 
£3 displacements of this work. 
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1 Introduction 

The problem of random nights has been appealing for many researchers in different scientific 
i-H fields. The original formulation is due to the statistician Karl Pearson, who, in a brief letter 

quoted in Nature, 1905, wrote: "A man starts from a point O and walks a yards in a straight 
line; he then turns through any angle whatever and walks another a yards in a second straight 
line. He repeats this process n times. I require the probability that after n of these stretches he 
rN is at distance between r and r + Sr from his starting point O." 

Pearson's aim was that of modelling the random migration of mosquitos invading cleared 
jungle regions, while Rayleigh in the same issue of Nature, observed that Pearson's problem is 
equivalent to the problem of the superposition of n sound vibrations with unit amplitude and 
arbitrary phase. The Pearson walk was generalized by Kluyver (1905) who considered steps 
with arbitrary but deterministic length. Successively, Rayleigh (1919) extended isotropic planar 
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random nights to the space R 3 , useful as a possible model of statistical mechanics of a diluted 
solution of polymeric chains. These random models also emerge in astronomy to describe the 
stellar dynamics as noted by Chandrasekhar (1943). The author also pointed out the link between 
random flights and diffusion processes. More recently, Stadje (1987) and Masoliver et al. (1993) 
dealt with a two-dimensional random walk moving with constant velocity and with directions 
uniformly distributed in [0, 2tt]. 

Over the years many papers, particularly in the physical literature, analyzed the properties 
of these random models, see Section 2 in Hughes (1995) and references therein. 

For the position X d (£), reached at time t > by the random flights in M. d ,d > 2, the 
conditional distribution 

P{X d (t) G d^\N d (t) = n} (1.1) 

has been the main object of investigation, where Afd{t),t > 0, is the number of changes of 
direction recorded up to time t. In the first part of the paper the number Afd(t), t > 0, of changes 
of direction up to time t is assumed to be a fixed number n. Franceschetti (2007) has obtained 



a condition for ( 1.1 ) to be uniformly distributed and this relates the dimension d of the space in 
which X d (t), t > 0, develops and the number n of changes of direction. A similar work has been 
carried out by Garcia-Pelayo (2008) and discussed by Le Caer (2010) who also considered random 
flights with Dirichlet distributed displacements. Orsingher and De Gregorio (2007) tackled the 
problem of random flights in higher spaces by dealing with uniformly hyperspherical distributions 
of the orientation of motion. Beghin and Orsingher (2010) considered a planar random motion 
where the deviations occur at odd-order Poisson events (or at even-Poisson events thus implicitly 
assuming that the displacements take a Dirichlet distribution). 

In the original Pearson's formulation of the problem of random flights the length of steps was 
deterministic and the probability distribution of X d (£), t > 0, lead to integrals which could not 
be explicitly worked out (see Watson, 1922, pag.421). For random flights in M. d with uniformly 



distributed steps it was possible to obtain the distribution (1.1) for d — 2,4, and for arbitrary 
values of n (and also the non-conditional distributions) . 

We remark that in all papers mentioned above the deviations are separated by exponentially 
distributed time lapses. This corresponds to assuming instants of changes of direction uniformly 
distributed under the condition that the number of Poisson events N(t) is fixed. This basic 
assumption permits us to obtain the explicit distribution of the position of the moving particle 
only in the spaces K 2 and K 4 . In R 2 the distribution was obtained by recursive arguments 
by Stadje (1987), by Masoliver et al. (1993). In R 4 the explicit distribution was obtained in 
Orsingher and De Gregorio (2007). This unlucky circumstance (for concrete purposes the space 
R 3 is clearly the most important one) is here overcome by the assumption that the intervals 
between successive changes of direction exhibit a Dirichlet distribution. 

By choosing a suitable basic parameter of the Dirichlet random displacements for each d, we 



are able here to obtain the explicit distribution ( 1.1 1 of X d (i), t > 0, for all n. Furthermore, the 
distribution has the universal isotropic form 

Mx d )^(c 2 t 2 -||xJ| 2 )" (1.2) 

where b depends on n and d, while A is the necessary normalizing factor (depending on t) . In our 
view this solves the classical Pearson's problem of random flights for all values of d in Euclidean 



spaces. We show below that functions of the form (1.2) are solution to some d-dimensional 
telegraph equation. 

By resorting to fractional Poisson processes we can randomize the distribution with respect to 
n, thus arriving at the non-conditional distributions of X rf (t), t > 0, for all d > 2. Furthermore, 
we can extract the distribution in R 2 and R known so far as particular cases of those derived 
here with Dirichlet distributed steps. 



We now describe the random flights analyzed in this paper. A random walker, starting from 
the origin of a frame of reference, moves in d-dimensional real space, with d > 2, at finite speed 
(denoted by c) according to the following rules. We assume that in the time interval [0,i], n 
changes of direction of motion are recorded. We suppose that the instants at which the random 
walker changes direction are < t\ < t 2 < ■ ■ ■ < t n < t, n > 1, and denote the length of 
time separating these instants by Tj = tj — tj-i, 1 < j ' < n + 1 with to = 0, t n+ i = t. Each 
displacement has an orientation defined in R d by the angles (9\,02,---,Od-2,^>) and we suppose 
that < 9j < n, < <j) <2n with joint law equal to 

g{0 1 , ....,$ d - 2 , 4>) = T ^ 1 l sin d - 2 B x sin^" 3 8 2 ■ ■ ■ sm9 d _ 2 . (1.3) 

2-7T2 

The law ( |1.3[ ) tells us that the direction is uniformly chosen on the hypersphere of M d with radius 
one. An important assumption is that the random vector (n, ...,T n ) (representing the length of 
the n displacements) possesses joint density equal to 

r((n + l)(d-l)) 1 "ft d _ 2 

Jl\Ti,...,T n )— (r(d-l)) n+1 t(»+l)(d-l)-l 11 ■? ' y ' 



where < tj < t — J2k=o T k> 1 — 3 — n > an d T n+ i = t — X)i=i T i- ^he distribution (1.4) 



is a rescaled Dirichlet distribution, with parameters id— 1, ...,d— 1), d > 2. The probability 



distribution (1.4) can be obtained as a marginal integral of the multiple uniform law as shown 
in Lachal et al. (2006). 

We will treat also the random flights with intermediate step length having joint distribution 

f( , r((n+i)(f-i)) 1 "+ 1 g _ 2 

h{n,...,r n )- (r( |_ 1))n+1 t(n+ i )( |_i)_i ll r j • t 1 - 5 ) 

where < Tj < t — X)fe=o T k i 1 — 3 — n i an d t„+i = t — 52,- = i T j > which is a Dirichlet distribution 
(suitably rescaled) with parameters (| — 1, ...., | — 1), where d > 3. 

The model treated here consists of the triple {6, t, J\f d (tj) of independent vectors where 9 = 



(6i, ...., 9 d -2, (f) is the orientation of displacements (with uniform law (1.3)), r = (ti,...,t„) 



represents the displacements and Af d {t) is the number of changes of orientation. In the models 



analyzed in Stadje (1987) and Orsingher and De Gregorio (2007), has law coinciding with ( 1.3 1, 
t is uniformly distributed and Af d (t) is a homogeneous Poisson process. 

Random processes with intcrtimes with non-uniform distribution have occasionally been con- 
sidered in the literature (see Di Crescenzo, 2002 and Pogorui and Rodriguez-Dagnino, 2005). 
Random motions in R c with n + 1 direction have been analyzed by Samoilenko (2001), Lachal 
(2006) and Lachal et al. (2006). 

Recalling that the motion develops at constant velocity c, the process described by the d- 
dimensional random flight X rf (t) = (X±(t), ...,X d (t)),t > 0, has components equal to 



n+i 

,-. dn a. . c.;^, a„ . ... dr. a . . . o;», j, . 



X d (t) — c y Tj sin 9ij sin 02j ■ ■ ■ sin 9 d -2j sin 

n+l 

X d -i(t) = c/ Tj sin 9\ j sin#2..j ■ ■ ■ sin0 d -2,j cos 




Figure 1: A sample path of a random flight in ]R 3 consisting of four displacements with angles 
(Oj,(f>j) defining each direction. 



(1.6) 



n.+ l 



X%{t) = c 2^ Tj sin 9ij cos #2, 
3=1 

ra+l 
Xi(t) = C 2^ Tj COS 0\ J. 



i=i 



If we consider a random flight in M. d with time intervals between successive changes of ori- 



entation distributed as (1.4), we will show that the the position of the moving point X^W = 
(Xi(t), ...,Xd(t)) at time t, has distribution 



Px rf (x d ,£;n) = 



r(=±i(d-i) + !)(c?v 



|x d ||2)f(^D-l 



r(f(d-l)) ^d/a^Cn+iXd-i)- 



(1.7) 



with <i > 2, ||x d || < ci and n > 1. For rt = the moving particle's position is uniformly dis- 
tributed on the surface of the ci-dimensional hypersphere H^ t with radius ct, and zero elsewhere. 



We remark that for d = 2, we can extract from (1.7| the planar distribution 

n 



Px,(x 2 ,t;n) 



27r(rf) r 



(cH 2 



|2\§-1 



|X 2 | < ci ' 



see formula (1.1) in Orsingher and De Gregorio (2007). We also observe that (1.7) has the 




Figure 2: The sample path of Figure 1 with its projection on M 2 (dotted line). 



remarkable structure (1.2) for b = %(d— 1) — 1 and the normalizing value equal to 



A = 



1 



r(=±i(d-i) + i) 



^d/-2^(n+l){d-\)-l r(f(d-l)) 



for all dimensions d > 2 and all numbers of changes of direction n. 



Furthermore, if we assume (1.5) as the joint distribution for the time length intervals (r 1; ..., t„), 
the position Y d (£) = (Yi(i), ..., Yrf(i)), t > 0, d > 3, of the random flight has probability density 
which reads 



PY d (y d ,£;n) = 



r((n+l)(|-l) + l)(c 2 i : 



2j-2 



2)n(f-l)-l 



r(n(|-l)) 7r d/2( rf )2(„+l)(f-l) 



(1.8) 



with d > 3, ||y || < ct. For d = 4, formula (|1.8[) simplihes into 



p-Y,{y,,t;n) 



n(n+ 1) 

7T 2 (C<) 2 "+ 2 



(cV-||y.|| 2 )' 



which coincides with the result (3.2) of Orsingher and De Gregorio (2007). 

The technical reason for which this important and simple result is possible is due to the semi- 



group property of the Bessel functions (2.7) (applicable with Dirichlet distribution (1.4)) and 



(2.9) (applicable in the case of Dirichlet distribution (1.5)). If we do not harmonize the order 



of the Dirichlet distribution with the dimension of space d, we get entangled in highly compli- 
cated formulae as Beghin and Orsingher (2010) showed. Throughout the paper the orientation 
(with distribution (1.3)) and the intertimes (with probability laws (1.4) and ( |1.5[ )) are assumed 
independent. 



We also observe that (1.8) corresponds to the uniform law if n = -jzv an< ^ tlh s implies that 



for d = 3 we need two changes of direction in order to obtain the uniform law, while for d = 4 



one change of orientation leads to the same result. In the previous case (1.7) the condition to 
have a uniform distribution is n = j 2 ^ and this means that for d — 2, n must be equal one in 
order to obtain this strange and unexpected result. 



For n = -^j the law (1.7) reduces to the form 



Px„ x d ,i 



d-1 



TT d / 2 {ct) d ~ 1 r(|) 



y (c 2 t 2 



(1.9) 



and this is valid for d = 2. For n = -^^ the second law ( 1.8 ) takes again the form ( 1.9 1 and this 
holds for d = 3. Therefore for these two cases the distribution displays a singular behavior near 
the surface of the sphere. In general, the bigger the number n of changes of direction the more 
concentrated around the starting point the distributions (1.7) and (1.8) are. This is because 



the sample paths coil up around the origin since they are subject to contradictory, fragmented 
displacements. 



From (1.7) and (1.8) we can extract the distributions of the random flights on all subspaces 
l , 1 < to < d, which preserve the structure (1.2). For to = 1, we are able to obtain the 



distribution of the projection of the random flight on the line by means of order statistics (as 
previously elaborated for telegraph processes, see De Gregorio et al., 2005, and some planar 
extensions with a finite number of directions, Leorato and Orsingher, 2004). The marginal 
distributions of (1.7) and (1.8) can be interpreted as the probability law of the random flight 
described by the shadow on the subspace E m , 1 < m < d, or, equivalently, as a motion where 
changes of direction imply also a random change of the velocity in the subsequent displacements. 
The projection of the distributions (|1.7[) and (1.8) on the line yields 



fj Cl {x u t;n) = 



r((n + i)(rf-i))(c 2 ^ 2 -x 2 )^( d - 1 )- 1 
(r(=±i(d-i))) 2 (2ci)(«+ 1 )( d - 1 )- 1 ; 



with d > 2,\xi\ < ct, and 



fy^yi^n) 



r(2(n + l)(f - 1) + 1) (cH 2 - y 2)( n +i)(#-i)- 



(r((n + l)(f-l) + |)) 2 



(2c*) 



2(n+l)(|-l) 



(1.10) 



(1.11) 



with d > 3, |yi| < ct. Particularly interesting is the case d = 3 in (1.10), because it yields the 
conditional distribution of the telegraph process T(t),t > 0, that is 



P{T(t) € dx\N{t) = 2n + l} 



„2\n 



(2n + l)!(c 2 t 2 -a: 2 ) 
(n!) 2 (2ct) 2n + 1 



\x\ < ct, 



where N(t), t > 0, is a homogeneous Poisson process. From ( |1.11[ ), for d = 3 and n = 2k— 2, k > 1, 
we arrive instead at 



P{T(t) G dy\N(t) = 2k} 



{2k)\ct (cH 2 -y 2 ) 



,,2\k-l 



\y\ < ct, 



fc!(fc-l)l (2ct) 2k 
which coincides with formula (2.18) of De Gregorio et al. (2005). For odd values of d we can 



derive from (1.10) the distributions of reinforced alternating processes on the line described in 



De Gregorio et al. (2005). 

The results ( |1.7[ ) and (1.8 1 depend on the dimension d of the space R d and on the number 
n of changes of direction. In order to obtain unconditional distributions (as in K 2 and K 4 in 
earlier work see Stadje, 1987 and Orsingher and De Gregorio, 2007), we here assume that the 
number of changes of direction is randomized and has the structure of a fractional Poisson process 
independent from the Dirichlet r.v.'s representing the step lengths and independent also from 



the r.v.'s representing the orientation of each displacement. Therefore, we average distribution 

owing distributioi 

P{N d (t)=n} = 



( 1.7) by means of the following distribution of a fractional Poisson process Afd(t), t > 0, 

1 (At)' 1 



E^A\t)T((^)n+iy 



(1.12) 



with A > 0, d > 2, n — 0,1,..., while for (KL8J), we must take the fractional Poisson process 
Md(t),t > 0, with probability distribution 



P{M d (t)=n} = 



1 



(xty 



E* hi (\t)T((l-l)n + iy 



with A > 0, d > 3, n = 0, 1, .... By combining ( 1.7|) and (|1.12| we obtain the probability law 
P{X d (i)edx d } Ai( c 2 i 2 -||x d |'^ 



./ -1 _-. E d— 1 rf-1 



At(c 2 t 2 -||x d || 2 )V 



n,=i^i 



n i {ct)^- 1 )- 



E d-i d (At) 



which, for d = 3, simplifies into 

P{X 3 (t) € rfx 3 } 

nLi dx 3 



x e A( c2 * 2 -i^ii 2 ) 

Tri cH 2 ^1.3 (At) : 



which is similar to the unconditional distribution in M 4 obtained in Orsingher and De Gregorio 
(2007). The functions appearing in the formulae above arc called two-parameter Mittag-Lefncr 
functions 

00 k 

^ = giW^)< *^>o,/?>o, 

and play a central role in fractional calculus. 

2 Exact probability distributions for a random flight in M. d 

We start our analysis by presenting the characteristic functions of the vector processes X d (4), t > 
0, and Y d (t),t > 0, defined as (1.6), when the number of displacements is fixed and equal to 
n + l,n > 1, and the changes of orientation are separated by random times ti,...,t„ with 



distribution (1.4) and (1.5) respectively. We denote by a d = (aj., ..., a d ) & K d while ||x d || = 



'X) 7 '=i x ? an d < a rf ,Xrf >= 2~Zj=i a j x i represent the Euclidean distance and the scalar product, 
respectively. 

Theorem 1. For the vector process X d (£) = (X\(t), ..., Xd(t)),t > 0, with intermediate time 
lengths with joint distribution (1.4 1, the characteristic function reads 



£ j e *<«d,X d (t)>j 



J ' ' -^J^-^MI^II), (2.1) 



(c*ll"dl 



;(<*-!)- 



where d > 2. For i/ie vector process Y d (£) = (Yi(t), ..., 1^(4)), £ > 0, to£/i intertime lengths having 
joint distribution (1.5 I, the characteristic function is 



*{—->} - 2W, t::!«:^ i)+i) w^(*faii). <-) 



MKII) 



with d > 3, where 



oo 

J„(a0 = 5>l) fc (f ) 



2k+u 



k=0 



k\r(k + v + i) 



x, v e e, 



is t/ie Bessel function. 



Proof. We show that under the assumption that (1.4) represents the joint distribution of the 
intervals n, ..., t„, the characteristic function of the position of the d-dimensional random motion 
X d (i) is equal to (2.1 ). We can write that 

rt—Ti 

'o Jo 



E 



, _ y^ n - 1 

^<a d ,x d (t)>\ = I dn I 1 dT2 . . . J J =' J ^ /i(t ^ ^ Tn )T n {a d ) 



where 
x n{a d ) 



(2.3) 






p7T pTT p7T p27T /*27T 

dO 1,1 ■ ■ ■ d9i, n+ i ■ ■ ■ dQd-2,1 ' " / d6d-2,n+i \ 
lo Jo Jo Jo Jo 

x I I < exp{icrj (ad sin 0\,j sin 8 2 j ■ ■ ■ sin dd—2,j sin 0j + a<i_i sin 0! j sin # 2 ,j ' ■ ■ sin &d-2,j cos </>j 



a 2 sin 6i tj cos 6» 2 j + a x cosfli,.,-)} d/2 sinfl^- 2 ■ ■ • sin6 d _ 2 ,j 



The multiple integral (2.3) is performed with respect to (n+l)(d — 2) angle variables ^j and 
n + 1 variables <j>j, with 1 < j < d — 2, 1 < j < n + 1, appearing in the orientation distribution 
( 1.3 ). The integral I n (a d ) has been worked out by Orsingher and De Gregorio (2007) as follows. 
Since 



1 

27 



2/r 



s ix(a cos <f+b sin 0)^ = J Q ( xy / a 2 + tf^ 



we observe that, after integrations with respect to <f>j,j = l,...,n+ 1, (2.3) becomes 

/■7T /*7T /*7T 

G$l,l ' ' ' / G?#l,n+1 ' ' ' / dQd-2,1 ■ ■ ■ d6d-2,n+l 

/o Jo Jo Jo 

n+1 r 

x 1 f < exp{icTj (ocd-% svnOi j sin6* 2 j • ■ ■ sva9d-2,j cos0d-2,j + ■ • • + a 2 sin^ij cos02,j + a i cos ^i,i)} 

r(d/2) 



x Jo (cTj sin Si j • • • sin9d-2,jJ(Xd + a d-i) rf / 2 _i sm ^ij 2 ' ' ' sin ^ 



smV d -2j 



We are able to perform all the (d — 2)(n + 1) integrations with respect to the angles 0i,j, 1 < 
i < d— 2,j = 1, ..., n + 1, by applying successively the formulas below 



"/ 2 , , fra v J v+ % (VaJT¥) 
(sinx) l/+ cos(6cosa;)J ! y(asina;)(ix = W — ; . 

(i V 2 (a 2 + b 2 ) 2 + 4 



(2.4) 



for Re v > — 1 (see Gradshteyn-Ryzhik, 1980, pag. 743, formula 6.688.(2)). The integration with 
respect to 0d-2,j, j = 1, ..., n + 1 yields 

r^ rf -2,i • • • l* Ma-2,M-i IT e ier i a *-*'* aB ^~'* a6 *->>''*» e *-*>i sin0 rf _ 2>j 

Jo Jo •_! 

Jo (erj sin 6»i j • ■ ■ sm6 d -2,j^a 2 d + a|_ x J 

g , CTj a d „ 2 si„ ei , J ...si„ 9d „ 3 , 3 co sed _ 2 , sin ^_ 2j J (cr, Sin^,, • • • S hl d _ 2 j^ + O^j) <%-2J 



n+l 



n 

7=1 

1+1 ( /.7T/2 

nW cos( 



CTjad-2 sin 6»i j • • • sm6d-3,3 cos6 d -2, 3 ) 



sia0d-2,jJo (cTj sin 6<i j • • • sm6d-2jJa% + ol\_ x ) dQ d -i 




^Ji/2 [CTj sin 0i ,j • • ■ sm6d-3,jyja d + a 2 d _ x + a|_ 2 

r/2" 



ct j sinWij • • • smC/ d _ 3 ,jWa^ + a 2 d _ x + a 2 d _ 2 



} . 



In the last step we applied formula (2.4) for 



v = 0, a = ctj sin 0i j • • sin 6 'd-3,j y a d + a d _ v b = cT. } a d -2 sinflij • • sin0d_ 3 j 
and also considered that 

/ sin(/3cosx)(sina;) ly+ J„{asvax)dx — 0. 
Jo 

The integration with respect to the variables 9d-3 i, ■•■, $d-3. n+i follows similarly by applying 



again (2.4) and yields 



/ d6d- 3 ,i ■■■ d6 d - 
Jo Jo 



n+l 

n„icTja d _ 3 sinSi^-- sin9 d „ 4 , : ,- cose d _ 3iJ - 

3=1 



sin! 



d-3,j 



x 2 



7T J V2 ( 



ct j sm f i,j ' ' sm Vd-3,j \ oq + a d _ x 



a 



(ctj sin 6*i j ■ ■ sm6d-3,jJa d + a 2 l _ 1 + a 2 



d-2 

r/2 



d-2 



n+l 



n(2 

i=i 



-\ 2 Ji ( CTj sin 6»i,j • • sin B d -Aj Jaj + a 2 d _ 1 + a 2 d _ 2 



l d-3 



ctj sin6»i,j • • sine/d-iWG! 2 , + a^_ a + a^_ 2 + aj_ 3 



By continuing in the same way, that is by applying successively formula (2.4) we obtain that 

(2.5) 



^,<> ^ r u n 7 



= 1 (CTj | lOrfl 



Therefore, the characteristic function becomes 

n+l ,t 



^| e *<sd^,(t)>| = M- 1 ! 



2)1 i Tl dT 2 ---\ dT n fl(n,...,T n 



^ "+} J< _i (cry | |a d | |) 
X M K-HaJI)^ 1 

, lr/ 'd\r i+1 r((n + i)(rf-i)) 1 



'"' ' ' '' (r(d-l))«+ 1 t(n+l)(d-l)-l 



/ 7f- a *i / r 2 d - 2 dr 2 • • • / r r f 2 (i - £ r^dr, 

Jo Jo Jo =1 

"+ 1 J 4 _ x (C7J | k,P 



X H toll^H)*- 1 (2 - 6) 

In order to work out this n-fold integral, the following result (see Gradshteyn-Ryzhik, 1980, pag. 
743, formula 6.581(3)) 

^(«- I )%(x)J,( a -# ^^ ) ^^ ) B ^iJ^.( B ), (2.7) 

o v27iT(/x + v + 1) ' 2 

with Re /i > — 2 and i?e ^ > — » , assumes a crucial role. Indeed, we apply recursively the formula 



(2.7) to calculate each integral with respect to the variable tj. In the first step we have therefore 









J|_ 1 (cr„||a d ||)J|_ 1 (c(i-^r ;7 -)||a d ||) = (cr„||a d || = y) 



1 /•c(*-E7= 1 1 T ; ,-)ll« d ll 

(c||a d n< 



/ dw 5 - 1 (c(t-E T i)Hsall-l') 5 

JO _•_, 



2d-3 

ra-l 

x Jd_ x (y) Jd_ x (c(t - Y^ T j)\hd\\ ~ V) 
i=i 

= . ,, „.-. , V- 2 jj (cft-yT ? )||a d ||) rf -^J d 3 (c(t- Vr 3 -)||a d l 
The second integral is given by 

( r(^i)) 2 i p-Y$3-r, 



2¥r(d-i)( c ||a d ||) M - 5 7o dr "^ 



:(cr n _i||a d ||)= 1 (c(* — J3 T i )||a d |]) d ' J r |_i(cr„_i||o cI ||)J (J _|(c(t- J2 T M&d\\) = (^nlk,!! = y) 

(r(^i)) 2 i /•c(*-z;?-ir,)ii a<J ii 



i-2 



27rr(d-l)(c||a d ||)3^7 dy 
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n-2 



n-2 



xy 



*-Hc(t-Y, T ^\\-v)^J$-i(v)J^$(<t- y ETj)\\<**\\-v) 



3 = 1 



i=i 



(r(^)) 



(V2^) 2 T(Ud-l))(c\\a d \ 



\3d-4 



n-2 n-2 



By considering formula (2.7|, we see that after (n — 1) integrations (with n > 2), the exponent 



of c(t — Ti)||a d || as well as the order of the Bessel function is given by the formula 



f-l)(n-l) + i(n-2) 



1 n , , , 1 
2=2^-^-2 



with /i = | — 1, ^ — (| — l) (n — 1) + |(n — 2). The exponent of c||a d || can be obtained by 
adding to dn — (n + 1) (where n here is the number of integrations) the number 2(f — 1) due to 
the adjustements necessary to apply (2.7). 
Then, the last integral becomes 



(r(W 



1 



)(d-l)(n+l)-2 



dr\ 



(2.8) 



(y2V)«-ir(f(d-l))(c||a d | 

x{cr l \\a d \\)^- l {c{t-T l )\\a d \\)^ d - l ^J i _ 1 (^i\K\\)J^ d ^ ) ^(c(t-r l )\\a d \\) = (ar 1 \\a d \\ = y) 



r 1 ))" 



c*l|a d ll 



dy 



(v^F)"- 1 r(§(d- i)) (c||a d ||)(<i-i)(«+i)-i 7 



(r(¥)) 



n+l 



(v^7F) n r(^i(n + i)) Mted 



1)(n+1) _ 1 J^i (d -i)-i(cf||a <i ||) 



ll) (d " 



Therefore, plugging the result (2.8) into the expression (2.6), and by observing that by means of 



the duplication formula we have that 



T\-\=^- dT{d - l) 



r(*=i) 



and 



r.-fc±a tf -i )+ i 



/7T2 1 



- ( n + i)( d -i) r((n+l)(d-l)) 

2 



r(fc+ii (d -i)) 



some simplifications lead to result (2.1) 



Under the assumption that the f%{r\, •••,T n ) is the density law for the intertimes Tj,j = 
1, ..., n+1, and by using arguments similar to those of the first part of the proof, the characteristic 
function of Y d (i),t > 0, assumes the following integral form 



E 



| e *<s (i ,Y d (t)>| 



d i 

22- x r 



n+l 



r((n + l)(f-l)) 



f 



(r(| - 1))™+! t (»+l)(f-l)- 



2-2 



dr\ 



t-n 



2—2 

t 2 2 dr 2 • • • 



t-E^ 1 



"A 1 J 4-i( 



CTi \\a„ 



f-2,, v^ ^-2, TT "v-\\^'i\\^d\\) 
r„ 2 (*->>■) = 2 dr„n 2 

j=i J= i (ctj I |Od| I) » 
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The first integral with respect to r n becomes 

rt-I2?=tTj J 4 _ 1 (cr„||a rf ||) J<f_i(c(*-X;™ =1 T.j)||a 



(c\\a d 



-4 / _ ll„. II „/-* V-« ,.MI„. II dT n = {y = CT n \\a d \ 



II 



cr„||a d | 



=(*-E?=i r j)llsdl 



n-1 



1 r (*-E^r J )|| Sd || j # _ l(y ) J # _ l(c ( t _ £™~ V^HodH -y) 



y 

in-1 



c(*-E?=lTj)lled!ll -2/ 



rf/y 



(c|KII) rf - 3 7 

1 2 J 2(|-l)( C (* _ E"=l 7j)ll«d 

= (c||« d ||) d - 3 |-l c(i-E-=iS-fcll 

where in the last step we have used the following formula (see Gradshteyn-Ryzhik, 1980, pag. 
678, formula 6.533.(2)) 



a JJx)J v {a — x) 



,, „ Ldx=[ l + l\J^ v {a) 



/i v J a 



Re n > 0, Re v > 0. 



/ x(a — x) 
The second integral provides us 

I 2 ft-^!^^ J r |-i(crn-i||Sdll) ^ 2 (f-i) ( c (* - Z)"=l TyOll^dl 



(2.9) 



vn-l 



(c||a d ||)t a -° f-Uo cTn-i||a d 



c (*-ELi7j-)l|odl 



i 2 /-^-Epi^oiisdii J^^y) ^2(|-i)( c (* — S"=^Tj)ll«dll -y) 



rfT n _i = {y = cT n \\a d 



^n-2 



(c||a d ||) 
1 



P- 4 f-l^u 



y 



c(*-Ei=i' r i)lladll -v 



dy 



n-1 



3 J3(|-l)( C (*-E"=l7-; 



jJW&d 



(c||a d ||)l d -M|-l) 2 ^-E^r.OII^II 



In the last integral, the exponent of c||a d || is equal to (n — 1 + 2)(| — 1) — 2 (n — 1, with 
n > 2, is the number of integrations performed). Therefore 



1 



J$_ 1 (cn.\\a d \\) J n (i-i)( c ( t - T i)\\&d\ 



(c||a < ,||)(»+ 1 )(#- 1 )- 2 (2- 1 ) n " 1 ^o CTi\\<*d\\ c(t-n)||a d || 

1 n /• ct II^IIJ|_ 1 (2/)4(f-i)( rf ll«dH-y) 



(c||a d ||)(" +1 )^-i)-i(f-l)™-iy y ct||a d ||-y 

1 U+l J (n+l)(f-l)( Ci H«dll) 

(clladlD^+^-D-id-l)" cfe|| 



-dm = {y = cn\\a a 



dy 



and then the result (2.2) follows immediately 



□ 



In the next Theorem we are able to invert the characteristic functions (2.1) and (2.2). 



Theorem 2. The probability laws ofX d (t).t > 0, and "Y d (t),t > 0, are respectively equal to 

(2.10) 



^fe^)^^- 1 ^'^-"^-''-' 



r(f(d-l)) ^/a^n+iW-i)- 
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d > 2, and 



px*<3U>*>*> = r(n(f-i)) 7rd / 2(d) 2(n + i)(f-i) • < 2 - n ) 



d > 3, mi/i ||x d || < ct, ||y || < ct and n > 1. 



Proof. By inverting the characteristic function (2.1 1, we are able to show that the density law 



of the process 'X. d (t),t > 0, is given by (2.10). Therefore, by passing to the hypcrspherical 
coordinates, we have that 

px d (x dl (;n) = -^ f e- i <^>E\e i< ^<^ t »}da 1 ---da d 



(27T) 

— ^ / /-Mp / <#!•••/ dOd-2 I d^sm^Ox-- -811x9^2 

( zn ) Jo Jo Jo Jo 



2tt 



x exp{— ip(xd sin 9\ ■ ■ ■ sin &d-i sin 

f^Hrfffd-iJ + i) 

(ctp)"" 



X2 Sin 6*1 COSC/2 + Xi COS V\ ) 



Jn±L {d -i)-l{ctp) 



2 ^i(d-l)-i / n + 1 js ,oo J|_ 1 ( / t.||x d ||)J a ± 1( ,,_ 1) _i(ctp) 

2 n+i (d -i)-i r (»±l (d _ 1)+ l) 

(2^)^/2 ( c £)^(<*-i)-§||x d ||#- 1 ./o 

1 r(^(d-l) + I) _ 

7r d/2 (ci) („ + i)(rf-i)-i r(=±i(d - i) - | + 1) l IUd 






+iCW_11_d±l 



|2N^±i(d-l)- 



In the first step above we have performed calculations similar to those leading to (2.5) and then 

f-271 



dffi ■ ■ ■ / <20 d _ 2 

o Jo Jo 

exp {— ip(xd sin #i • • • sin Bd-i sin 
rf J"d_ 1 (p||x d ||) 



sin d 2 6>i • • • sin 6> d _ 2 



(2.12) 



X2 sin 0i cos #2+^1 cos Q\ ) } 



(2^ 



(p||x d ||) 4 



while in the last step we have used the formula (see Gradshteyn-Ryzhik, 1980, pag. 692, formula 
6.575.(1) with a correction in the bounds of p and v) 

(a 2 - B 2 Y-^B^ 

J u+1 {ax)JJpx)x^ v dx = { ^ - -, 

^ »+i\ ) mvm > 2 v ~i J -a v + 1 T{v - p + 1) 

a>P, Re{v + l) > Re{p) > 0, for v= ^(d- 1) - |, p = f - 1, a = ct and £ = ||xj|. 
Analogously, for the random flight Y<j(t), i > 0, we have that 

, x , 2("+ 1 )(i- 1 )r((n+l)(|-l) + l) f"> d _ 1 -/|_ 1 (p||x d ||)^i)(f-i)(ctp), 

PY rf (y,,i;n) = zrT ^ — ^ / p -^ — iisd - — w „ 2 , lWd 1N dp 



-d^d' 



(2tt)* 7o (pllx.lDf- 1 (ctp)^ 1 )^- 1 ) 

2 (»+i)(#-i)r((n + i)(f-i) + i) r 

(2 7 r)i(ci)(™+ 1 )(5-i)||x d ||i- 1 7o 



p |-("H-l)(#-l)j | _ l(/0 || Xd ||)J (n+1)(| _ 1)(cip)d/0 
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r((n + l)(f-l) + l) .^^a.n n^cf-D-i 



Kd/2 ( ci )2(n+l)(f-l) r(n( | _ X)) 



n 



Remark 2.1. FKe can also check results (2.10) and (2.11 1 by evaluating their Fourier transforms, 
thus showing that we reobtain results (2.1 1 and (2.2 1 respectively. For T-L^ t = {x d : ||x d || < ct}, 
we have that 



/ e i<s *' ;s * > px d (Xd,t;n)dxi--dxd 

r(2±i(d_i) + |) 



(2.13) 



7r d / 2 r(f (d- i))(ct)(™+ 1 )( d - 1 )- 1 y wi 
r(2±i(d-i) + |) 



s i<a <i ,x d (t)>( c 2 t 2 _ ||^|| 2) f (d-l)-l dxi ... dxd 



7r d / a (rf)("+ 1 )(« I - 1 )- 1 r(f (d - 1)) Jo 

x exp{zp(cid sin 0i • • • sin dd-i sin </>+■■■ + 

x (c 2 i 2 - p 2 )t (^-D- 1 sin d - 2 0x ■ • • sin0 d _ 2 



/ _1 dp / d»! 

o 



d6 d -2 

^0 



2 d/2 r (s±i(d-i) + i) 

(Ct)(n+l)(d-l)-lr(f(d-l))7 



•> - ■' - , pd - 1{c 2 t 2 _ p 2 )W -l)-l J i-^P\^\\) 



a.2 sin 9i cos 62 + a± cos 6*i ) 

2*11. . 
TZT d P> 



(p\h d \\r- 



where in the last step we have used the result (2.12). 
Now, we work out the previous integral 



Mla s ||) ! 



(2.14) 



(-1) A 



V^ (-1) ([[grill) / ,x(»+l)(ri-l)+2fc-l /" l+fc-l/j 

^A:!r(fc + f) 2 2fe +l W 7o 

r ^ (d _ 1) \g(-i) fc MKH) 2fc 



(Ci)(n+l)(d-l)-l 



-y) 



t(d-i)-i 



rf;y 



fc! 2 2fc r(fc+2±i(d-l) + §) 



function (2.1 ) 



Therefore, plugging in the expression (2.14) in (2.13), we immediately obtain the characteristic 



Similar calculations hold for the characteristic function of the distribution (2.11) and then it 



is not hard to obtain the result (2.2) 



Remark 2.2. It is not hard to show that px. (x rf ,i;n) and py (y,,t;n) integrate to 1. Indeed 



~ d 

\ p^Md^ n )\\ dx 

Jn% j=1 



2lW(d-l) + i) 



ct 



r(f)(ci)("+ 1 )( d - 1 )- 1 r(f (d- 1)) Jo 

(p = ctVv) 

r(^(d-i) + |) f 1 M t(d _ 1M 

r(f)r(f(d-i)) Jo v ' 



P d - l {ch 2 - p 2 )t( d -D- 1 d p 



dj/ = l. 



For the density law py (y,,t;n) similar calculations hold. 
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ofX. d (t),t > 0, and Y d (t),t > 0, respectively 



Remark 2.3. From (2.10) (for d — 2) and (2.11) (for d = A), we can extract the distribution 

Px 2 (x 2 ,i;n) = —^—{cH 2 - feint- 1 (2.15) 



2ir{ct) r 



"0 + 1) / c a t2 _M ,|2sn-l 

-i'-^-I 7r 2(' c ^2n+2 l - IIJL4IM 



PY.(y,,*;n) 



(2.16) 



which have been obtained in Orsingher and De Gregorio (2007) under the assumption of dis- 
placements separated by intervals with uniform joint distribution 



/(ti,...,t„) 



t" 



where < Tj < t — J2k=o Tk > •"■ — J — n ' T "+i = ^~ S?=i T i- ^ n ^is case an homogenous Poisson 
process governs the change of orientation of the steps ctj, j — 1, ...,n. 



Remark 2.4. From (2.10) and (2.11), we are able to derive the exact distribution of a random 
flight moving in IR 3 , i.e. 



Px,(x3,*;n) 



r(n + |)( c 2t 2 -||x 3 || 2 y 



T(n) 



TT 2 (ct) 



2n+l 



r(^i + i)(c 2 i 2 -||y 3 || 2 )t-i 



^3 ^3 ' 



n > 1, I1X3II < ct 



n> l,||yj| <ct. 



r(f) ^f(ci)«+i 

The following relationship between the distributions (2.17) and ( |2.18 1 emerges 

Px, fej f ; "-) = PX, (x 3 , i; 2ra), 



(2.17) 
(2.18) 

(2.19) 



which tells us that a random flight developing in M 3 according to the Dirichlet law /i(ti, ...,t„) 
has the same density of a three-dimensional random flight where the steps have joint distribution 
/ 2 (ti, ...,t 2 „). 



For n — 1 and n — 2, the distributions (2.17) and (2.18) provide us the uniform distribution 
inside the sphere *M? ct (as emerges from Table\Sty, that is 



Px,(x 3 ,£;l) =pY,(y,,i;2 



3 w_3' 



r(f) 

7T2 (ci) 3 



Furthermore, 



observe that 



px,(r 3 ,fA) = ^ w 



:H 2 



which corresponds to (4-la) in Orsingher and De Gregorio (2007) for n = 0. The random 
traveller for one change of direction is more likely to be near the sphere surface cW 3 f while for 
two changes of direction his position is uniformly distributed inside the sphere. This has been 
commented by Franceschetti (2007) for a planar random motion. 

Remark 2.5. It is interesting to note that originally the problem of the random flights has been 
tackled by considering the length of the steps ctj constant and equal to A. Then, in this case the 
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following distribution emerges (compare with Watson, 1922, pag.^21, with suitable adjustments 
of the parameters) 



Px d (*d,t;n) 



2 (!- 1 >- 1 7 r-«[r(f)] n+1 



(A||x„| 



2)5" 



/ j \ / \ n-\-i 

p-(^- 1 )" +1 j | _ 1 (p|tall)(j|_ 1 (pA)) dp 



The above integral can not be worked out. Therefore, as suggested by Theorem [J| it is crucial 
to randomize the length of the steps in order to provide a general solution of the problem of the 
random flights for the real space having dimension d. 

The distance from the origin of the position reached by the d-dimensional random flights 
X d (t),t > 0, after n + 1 steps, that is Rd(t) = ||X d (i)||,£ > 0, has the following distribution 
function 



P{R d (t) < r} 



l! "' l «' / 1 ) + l) [" p d-l {c 2 t 2 _ p 2^{d-X)-X dp 



( ci )(«+i)(rf-i)-i r(f)r(f(d-i)) J 

and therefore the density law of Rd{t) becomes 

t)d (r t-n) = 2 2 - — - - 

PRAr,hn) r( | )r( s( d _i)) (ct) (n + i)(d-i)-i 

with < r < ci. Analogously, for L d (t) = ||Y d (t)||, t > 0, we obtain that 

r((n + i)(f - 1) + 1) i^jcH* / 2 )"(i-D-i 
PLdl '' nj r(|)r(n(| - 1)) (ci) 2(n + i) ( |-i) 



(2.20) 



(2.21) 



with < Z < ct. The behavior of densities ( |2.20| ) and (2.21 ) is outlined in Figure [3] for d = 3 and 
different values of n. 

We present now the expression of moments of the radial processes Rd(t),Ld(t), t > 0. 

Theorem 3. For p > 1 and n > 1, we have the following general results 

Y(P±i)r(n+±(d- I) + I) 



and 



r(|)r(^ + f(d-i)) 



^^y- ^^^t-^t^ W <>* (2.23) 



r(f)r(*±4 + „(f-i)) 



Proof. In view of (2.20), we obtain that 



£{i? d (i)} P 



2 r(2±i(d-i) + |) 



r(|)(ct)(™+ 1 )( d - 1 )- 1 r(f (d - 1)) 7 



^H-d-1^2 _ r 2 )f (d-l)-l dr 



r(f)r(|(d-i)) jo 
r(^)r(^±i(d-i) + i) 



''(i-w) 



Kd-1)-1 



rf.y 



r(|)r(^ + f(d-i)) 

By similar steps we arrive at (J2.23 ) . 



{cty. 



□ 



1G 



n= i / 

n=2 / 

n=3 / 

n=4 / 

i / / y '• 

i : / * \ 

i : / / > ■. 

/ •' / / v '. 

/ ■' - / v \ i 

i ■' ■ / x '. ', 

/ / / v \ i 

i .■ / / \ . . 

'•■ -' -/ x • ', 

<•'-' -^ x \ i 



0.0 0.2 



0.6 0.8 1.0 




Figure 3: The behavior of the densities of R%{t) and L 3 (t), with c = 1 and t = 1, for n = 1, 2, 3, 4. 



Remark 2.6. From (2.22) and (2.23), we can extract the following results 

d 



E{R d {t)} 
for d > 2, n > 1, and 



iW; 



r(|) n(d-l)+d 



"' ' ct, S{i? d (t)} 2 



n(d- 1) + d 



(ci) 5 



E{L d (t)} 



r(%i) 



d-1 



T(f) (n + l)(d-2) + l 
/or d > 3, n > 1. M^e observe that 



E{R d (t)} 



1 



1 



<*, £{£.*(*)} 



E{R d {t)} 2 



i(d-2) +d 



(ci) 2 



d-l+^r' E{L d {t)Y 



1 



n+1 



d-1 



(2.24) 



(2.25) 



E{L d (t)} 
and this shows that 

E{R d (t)} < E{L d (t)}, E{R d {t)} 2 < E{L d (t)} 2 
for all d > 3. 

Theorem 4. The projection of the processes X d (t),i > 0, and Y d (£),i > 0, onto a lower space 
of dimension m, leads to the following marginal distributions 



fi(* m >t;n) = 



T(^(d 1) + |) (c^ - ll^ip)^^^-!)-^ 

r( ( " +1 ^ (d i) i 1-m ) t^ (ct)( n+1 )( d - 1 )- 1 
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(2.26) 



/■£ (y ,t;n) 



r((n + i)(f - i) + i) (c 2 * 2 - \\y m \\ 2 Y n+1 ^-v- ! 
r((n + l)(f - 1) + 1 - f ) Tr^Cct) 2 ^ 1 )^- 1 ) 



(2.27) 



to£/i ||x m || < ct, ||yl| < ct and 1 < m < d. For m = d t/ie densities (2.26) and (2.27) coincide 
with (ElO| anrf pTfTl). 



Proof. We start by observing that the projection of the random process X d (i),i > 0, onto the 
space R m , represents a random flight with m components having density law given by 



/x (Xm)*; n )= / , dx m+1 - 



-\A 2 * 2 -I|x m ll 2 



V C2 * 2 -H^-2H 2 /■V /c2t2 -|l^-lH 2 

dx d _i / Px (x d ,t;n)dx d . 

- V ' c2 * 2 -ll^-2ll 2 ^-V^ 2 * 2 -!!^-!!! 2 



Then 



2 



.Px d (x d ,<;^)da;d 

r(s±i(d-i) + |) /•V c2 * 2 -n^-iii i 



^^(^(n+ixd-ij-i r(|(d-i)) 



(cV-lfell 2 )*^ 1 )- 1 ^ 






^(^(n+lXd-lJ-l r(f(d-l)) 

i r(s±i(d-i) + |)V5F, 2+2 



5(f- W )t( d - 1 )- 1 d W 



7r d/2 (ct )(n+l)(d-l)-l r(f(d-l) + i) (C< ll^-lH ) 

For the integral with respect to Xd, we have that 



2sa ((J _i)_i 



r(^±i(d-i) + i)V7r /•\/ c2 * 2 -n^-2ii 



7 r d / 2 (ct)("+ 1 )( d - 1 )- 1 r(f(d-i) + |) 



v - iisd-iir)*^ 1 *-*^--! 



(z d _i = ^^c?t 2 - ||x d _ 2 || 2 ) 

1 r(^(rf - f) + |)0r 2 2 _ j|2Nf(d-i) r 1 i (1 _ )t(d -i)-i rf 

r( ^ (rf - 1) + ^( C 2 t 2 -||x^ 2 || 2 )t(-). 



7r d/2( ci )(„+i)(rf-i)-i r (s(d _ i) + 1) 

Therefore calculating the successive integrals in this way, we obtain that 

r(fi±l(d- 1) + i)7T^ (c 2 t 2 - ||x m || 2 )f (d-D-i+ s 



/x_(Xm.*;») 



r(f(d-l) + ^f 2 ) 7r d / 2 (c*)C n+1 K<i-i)-i 



with ||x m || < ci, and by simple manipulations the expression ( |2.26[ ) emerges. 

By using the same approach we derive the result (2.27) concerning the projection of the 
random nights Y d (t), t > 0. □ 

Remark 2.7. The functions of the form 

q{* d ,t) = {cH 2 -\\?L d \\T> ||x d ||<rf, 
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appearing in all densities (2.10), (2.11), (2.26) and (2.27), are solutions to the telegraph-type 
eguations 

d 2 Q 2Y^^9 2m—l-\-ddq 



dt 2 



^dx 2 



t 



dt' 



(2.28) 



In particular, for m 
equation 



-^2^ equation (2.28) simplifies and becomes the d-dimensional wave 



d 2 q 

dt 2 



f-fflx?' 



Remark 2.8. From (2.26), /or d — 2 and m — 1, we obtain the marginal density of the con- 
ditional distribution of a planar random flight with uniformly distributed switching times, that 

is 



f 2 Xi {x l ,t;n) 



r(f)r(f + i) /2 



2.2 



5 I < rt 



X, 



n > 1, |xi| < ct. 



27iT(n) 

Furthermore, by setting d = 4 and m = 1,2,3 m i/ie densities ( 2.26[ ) and (2.27), we derive the 
probability distributions of the projections of a four-dimensional random flight, respectively X 4 (£) 
and Y 4 (i), onto i/ie lower spaces (see Table\w. Analogously, Table\w summarizes the marginal 
density laws of the three-dimensional random flights onto the spaces R 2 and R. 

From Table\l\emerges that for n = 2r+l, r = 0, 1, ... i/ie distribution of a planar random flight 
(with intertimes having uniform law) px (X;2'^ n ) coincides with / x (x 2 ,£;n). Furthermore, the 
distribution pY (y . , £; n) coincides with the result obtained in Orsingher and De Gregorio (2007), 
formula (3.2). Therefore, the marginal densities o/Y 4 (t),t > 0, of the Table below coincide with 
the probability distributions (4-la), (4- lb) and (4-lc) in Orsingher and De Gregorio (2007). 



d = 4 
m = 3 
m = 2 
m = 1 



X 4 W 



(r(fn+2)) 2 2 3 "+ 3 , o 2 ,, i.aNfn-i 

i 2 (d)3"+2r(3n+4)i CI IIS3II ; 2 



jrt+1 
7r(ci)3"+ 2 



( C 2 i 2 -||x 2 || 2 )f 



r(f »+2)r(f n+l)2 3 " +1 , o 2 2\fn+A 

77(ct) 3 "+ 2 r(3n+2) ^ C C X l/" 



Y4W 



r(«+2)r(»)2 2 "- 1 / 2f 2 || v ||2>,n-J 

2nV\ c l -iiy 3 n J • 



7r 2 (ct) 2 "+ 2 r(2n) 



3 1 
2\ra 



> n Z L > „ fr- 2 r - llv ll 2 l 
7r(ct )2„+2 l^c 1 ||y 2 IM 



r(n+2)r(n+l)2 



2+2 _ „,2\n+A 



(TH-jji jn+l^ / 2^2 _ ,,2-1 

7r(ct) 2 "+ 2 r(2n+2) \ C b Hi) 



Table 1: The density laws of the processes representing the projections onto the lower spaces of 
X 4 (t),*> 0, and Y 4 (£),£> 0. 



We are able to provide an alternative stochastic representation of the motion on the real 
line representing the projection of the d-dimensional random flight X rf (t),i > 0, as follows. We 
consider a random walker moving, with velocity c > 0, forward for a time £T„, where T„ is a 
Beta( r!: Y^(d — 1), ^^(d — 1)) random variable, and backward for the remaining time t(l — T~). 
In other words, the position X(t),t > 0, reached by the random traveller at time t is equal to 



X(t) = ct [T* - (1 - T„ d )] 



(2.29) 



Therefore, the density law of X(t),t > 0, becomes 
d 



p x (x,t;n) 



dx 



P 



{X(t) < x] 



= —plr d < 2ct + x 

dx I " 2ct 
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d = 'S 
m = 2 
m = 1 



X3W 



Ya(t) 



n+J 



("+!) 



|2\^ 



^^(A 2 - ||x 2 |n«-i ^^r(c^ - ||y 2 |P)^ - 

r (»+f) (,,2,2 _ ^,2\n r (f + f) I J2J.2 _ ,,2\ 3 

^(dj^+irfn+l) ^ L l X l-' VS(o*) n+lr (? + 1 ) ^ ^' 



Table 2: The density laws of the processes representing the projections onto the lower spaces of 
X 3 (i),t> 0, and Yg(t),t > 0. 



r((n + l)(d-l)) d 
(T(n±l(d-lWdxJ 



w 



^^-^(i-wy-^^-^dw 



(2ci)("+ 1 )< d - 1 )- 1 (T(a±i(d-1))) 



r,( "+ li( ' / - 1 )) (c 2 i 2_ a; 2 ) a±l(d-l)-l 



(2.30) 



and coincides with the distribution (2.26) with m — \. We observe that from (2.30) for (i = 3 
and m — 1, wc have that 



/iifcijtsn) 



P{T(t) e dxi|JV(i) = 2n+l} (2n + l)! 



1 



2\n 



dxi 



(n!) 2 (2ci) 2 ™+ 



l(c^-^) 



where the duplication formula for Gamma functions has been applied and T(t), t > 0, represents 
the standard telegraph process. The above distribution has been obtained by De Gregorio et al. 
(2005) (see formula (2.17)) by applying the order statistics. 



We can represent the one-dimensional motion underlying distribution (2.30) as a sequence of 



alternating forward and backward displacements at speed c, where the change of direction occurs 
at Poisson times tj with Tj = tj — tj—i, j = 1, ...,n. Since the intertimes tj, j — 1, ...,n possess 
uniform joint distribution, in force of exchangeability we can rearrange the displacements and 
put together forward steps and backward ones. In our case, we have for d = 1r + 1, r = 0, 1, ..., 
{n + l)r forward displacements and an equal number of backward ones. The instant T?^\ 
indicates the time where the last forward displacement occurs in the rearranged sequence and 
corresponds to the (n + l)r-th order statistics from a uniform distribution in (0, £). This leads 
to the representation (2.29). 



The d-dimensional random flight producing in R the planar motion described, for example, 



in Stadje (1987), with distribution (2.15), must satisfy the following relationship 



r = ( n + l)(d- 1) -1 

among d, the number r of changes of direction n in M. d and the number of changes of direction 
of the probabilistically equivalent motion in R 2 (in ( |2.26 ) write (n + l)(d—l) = r + 1 and m = 2 
so that (2.15) emerges). 



Analogously, for the projection onto R 1 of Y d (£),£ > 0, a similar representation holds true. 
Indeed, in this case we can write 



Y(t) = ct ^ (1 - S^] 



(2.31) 



where S d is a Beta((n + 1)(| — 1) + ~, (n + 1)(| — 1) + |) random variable. Then, we obtain 
that 



PY-(y,t',n) = 



dy 



p { f v<y} = Ty p { s *< 



2ct + y 
2ct 



(2.32) 
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r(2(n + 1)(| - 1) + 1) (c 2 t 2 - y2) («+i)( l-i) -I 



(r((n+l)(f-l) + i))) 2 



(2ct) 



2(n+l)(f-l) 



which coincides with the distribution (2.27) with m = 1 



Remark 2.9. The density f^ (x m ,i;n) becomes uniform on the hypershere "H™ when d 



"~*~4?i i 1 < w < d. Analogously, the distribution (2.27) is uniform on the hypershere H™ if 



d = 2 "+"+ 2 i < m < ,f Admissible combinations of m and n are those for which the dimension 
d is an integer number (see TableUM. Some simulations presented in Le Caer (2010) give graphic 
evidence of the uniform law for n — 2 and d = 2. 



*«*(*) 



X d (t) 



n = 2, d = m = 2 
n = m = 1, d = 2 
n = 1, d = m = 3 



n=l, m = 2, d = 3 

n = 2, d = m = 2> 
n = l. d = m = A 



Table 3: The values of n, d and m leading to the uniform distribution. 



Remark 2.10. We expect that the density tends to infinity near the surface of the hypersphere 
H™ for a small number n of deviations and this is confirmed by the following inequalities n < 
— 7_ p , for Xrf(£),i > 0, and n < m 7_ t > f or Y d (t),i > 0. Since m < d and n is integer, we 
obtain that the first inequality holds for n — 1 and d = m = 2, and the second inequality is valid 
for n — 1 and d = m = 3. In all the remaining cases the distributions take a bell-shaped structure 
because the larger is the number of changes of orientation, the shorter become the displacements 
and the closer to the origin is the moving particle. 



3 Unconditional probability distributions 

In order to obtain unconditional densities for X d (t), t > 0, we randomize the number of deviations 
Afd(t), at time t > 0, by assuming that it possesses the distribution of a fractional Poisson process 
(consult on this point Beghin and Orsingher, 2009). In this context by fractional Poisson process 
we mean a process with distribution 



P{N d {t) = n} = 



1 



(At)™ 



E^ A {\t)T{{^)n+iy 



d> 2, n = 0,l, 



(3-1) 



where E a: p(x) = J2T=o TTak+B) > ^ G 1^7 a : /? > 0, is the generalized Mittag-LefHer function. The 
generating function of the probabilities is 



GNa {u, t) 



Ed-i d(Xtu) 

2 ' 2 

E d-1 d (At) 



\u\ < 1. 



Since 



—E v p(ax) = - [£ , „.^+/j-i(a.T) + (1 - f3)E vv+ p{ax)] 
dx v 



(3.2) 
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we obtain that 



E{N d {t)} = j-G^(u,t)|„=i 

au 

2At 



(d- l)E^i d (Xt) 
The above result can also be obtained directly by using 



E^^X^+^l-fjE^^Xt) 



(3.3) 



E{N d (t)} = 



1 



d- lEd^i d(Xt) 



rffd-1 



(\ty 



^((^ijn+l) 



d-1 



d d 

- -1 hi 

2 2 



and by performing some straightforward calculations. If d = 2 the mean value of Afd(~t) becomes 



E{M 2 (t)} 



2At 



%i(At) 
while for d = 3, we get that 



Ei i (At) = (from (3.2)) 



£{A/- 3 (t)} = Ai- 



i^dK i(Af) = 4 l0gi M At) 



At ^i,| (At) 



2 ^i,s (At)' 



(3.4) 



Result (3.4) shows that the fractional Poisson process with distribution (3.1 1 has a mean number 
of events growing more slowly than the classical Poisson one. 

Analogously, for Y d (t), t > 0, we represent the random number of deviations by means of the 
process A4 d (t),t > 0, having probability distribution 



P{M d (t)=n} = 



1 



(Xt) n 



^_i,$(At)r((!-l)n+f)' 

The generating function of the probabilities in this case reads 

Ed_ 1 d [Xtu] 



d > 3, n = 0, 1, 



(3.5) 



GM d (u,t) 



E d_-x d \Xt) 



u < 1. 



From relationship (3.2) emerges that 

— G Md ( u ^)\u=i 

2At 

(d-2)E i _ 1 AXt) 



E{M d {t)} 



^4-i, d - 2 (At)+(l-2)^#-w-i(At) 



which for d = 4 becomes 



E{Mi{t)} = Xt[l 



^1,3 (At) 



£i,2(At), 

For Af d (t) = the particle reaches the surface of the hypersphere with probability 

1 1 



P{X d (t) e m%} = 
while if Ai d (t) = 0, one has that 

p{Y d (t) g m%} = 



Ed^d(Xt)T(l) 



Ed.AXt)^) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 
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We assume that the processes Afd(t) and A4d(t),t > 0, are independent from the Dirichlet 
distributed displacements and the angle orientations. Now, we provide the main results of this 
Section. 

Theorem 5. // the number of deviations is represented by a fractional Poisson process with 



distribution (3.1 1 then the absolutely continuous component of the probability distributions of 



X d (t),t > 0, is equal to 



P{X d («)edxJ Ai(c 2 t 2 -||x d ll 2 ) 



_-. E d-i d-i 



At(c 2 * 2 -||x d || 2 )^r 

(ct)< d -!) 



Y\Udx 3 ** (<*)**- 



E d-i d (At) 



(3.10) 



where d>2, ||x d || < ct, while if the number of the deviations is given by (3.5) , the distribution 
o/Y d (t),t > 0, reads 



P{Y d (t)edy d }_ \t(cH*-\\y d 



|2\f-2 Ed_ 1 d_ 1 



(ct)( d - 2 ) 



UUdyj ^ ( ci ) 4(f ~ 1} 

where d>3, ||y || < ct. 

Proof. For the random flight X d (t),t > 0, we have that 

g{=# -. d ~ d} = E^tmlfW) = n} 

[[j—i<iXj n=1 

1 1 



E d i rf (^) 



(3.11) 



^ (At)" (c^-HxJI^f^- 1 )- 1 

2 ' 2 n— 1 z v y 



1 1 vA (At)" +1 

2 ' 2 n— v ^ v 



"-' (c^-UxJI 2 ) 5 * 1 ^- 1 )- 1 



2 ' 2 

1 At 



)) ( ci )(«+2)(d-l)-l 



c 2 * 2 -||x rf || 2 ) 



d — 1 _i OO 



7r# Ea=i d(Xt) (ct) 2 ( d " 



1)-1 Z-r r fn±lM_! 



2 ' 2 



(W-1)) 



.> .i . iiO\ d—1 i Hj d — 1 d— 1 



At (c 2 t 2 -||x d I H 
tt! (ct) 2 ^- 1 )- 1 



n=0 



At(c 2 t 2 -||x d || 2 )^ 

(Ct)^" 1 ) 



E d-i d (At) 



Similarly, for the random flight Y d (t), t > 0, we are able to derive the density (3.11 ) as follows 

P{X d (t)edy d } ~ , ,......, , 

nd - ~ d =}_^PY d (y d ,t;n)P{M d (t)=n} 

U j= i d Vj 



n=l 
1 1 



E 



(At)" (c 2 t 2 -HyJI 2 )"^- 1 )- 1 



n l Ed^d(Xt) ^ r(n(| - 1)) ( c< )2(«+i)(f-i) 
1 1 ~ (Xt) n + 1 (c^-HyJI 2 )^ 1 )^" 1 )- 1 

_4 E, . ,(\+\ 2.^ 



ni Ed_ h d(Xt) ^ r((n + l)(f - 1)) (ci) 2(n+2)(f-i) 



At (c 2 t 2 - | |y, 112 ^- 2 - 






^Ed^diXt) (ct)4(|-D ^r((n +1)(|-1)) 



At(c 2 t 2 -||yJ| 2 )^ 

( ci )(d-2) 
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Xt{c 2 t 2 -\\y d \\ 2 ) 



2^4-2 #4-1,4-1 



(ct)< d - 2 ) 



7T 2 (ct) 



4(|-1) 



Ed_ 1 d(Xt) 



D 



Remark 3.1. Ji zs noi /iard to verify that 

f p{x d (t) e <***} = 1 - £{x d (t) g aw^} = 1 

and 



E^AXt)T(i) 



[ P{Y d (t) G dyj = 1 - P{Y d (t) G 0W* } = 1 - - * -1 



Remark 3.2. We examine some particular cases where the densities of TheoremU)\take interest- 
ing forms. Indeed, by observing that En(x) = e x and Ei 2 (x) = e , we have the distributions 
summarized in the following Table. 







d = 2 


d = 3 




A E 






x*W 


1 i(^V c2 * 2 -ifeii 2 ) 

2 ' 2 V v ' 


X e ^< c2t2 -ll^3ll 2 ) 


ttcEi , (At) 

2 |X 


V c2 * 2 -ifeii 2 

d = 3 


d = 4 




AcH E \,l^) 






X d (t) 


^i,iUV c2 < 2 -iiy 3 ii 2 ) 
V c2 * 2 -iiy 3 n 2 


A 2 e ^(<= 2 ' 2 -Hy 4 ii 2 ) 


w 2 c 4 t 2 e A '-l 



/orm ( 1.2 I 



^4// i/ie distributions (3.10), (3.11) ond i/ie special cases in the above Table have the isotropic 



If we suppose that the changes of direction are governed by an homogeneous Poisson process 
the intervals Ti,...,t„, are uniformly distributed on [0,i\. In this case the absolutely continuous 
component of the unconditional distribution of a planar random flight X 2 (t),£ > 0, is given by 
(see Stadje, 1987) 



Px,fe,£) 



27TC A / C 2 t 2_||x 2 ||2' 



|x 2 || < ct. 



(3.12) 



Therefore, by comparing (3.12) with the corresponding distribution in the above Table, we observe 



that in our context, e xt /2 ande^y° 2t2 H- 2 " 2 are replaced by E 1 ^(Xt) and Ex x [\\/c 2 t 2 - ||x 2 l| 2 V 
respectively. 

Furthermore, the absolutely continuous part of the distribution of a four- dimensional random 
flight with Poissonian switching times (see formula (3. 7) in Orsingher and De Gregorio, 2007) 
is equal to 

X b-\\tr II 2 I - A . 

|y, 



% 4 (y 4 >')- cW 



e -^iiyJ 



2+^-(c 2 t 2 >- " 2 



c 2 V 



which has to be compared with 



x 2 e ;M c2 ' 2 -iiy 4 ii 2 ) 

PY 4 (y 4 ,0 - ^^2" ^xT^—i 
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We indicate with Xf n (t) = {X 1 (t), ...,X m (t)),t > 0, and with Y^(t) = (Yi(t), ...,F m (t)),t > 
0, the random processes emerging from the projection onto R m of X d (t) and Y d (t), respectively. 
In the next Theorem we give the unconditional distribution of X m (t),i > 0, and Y m (t),i > 0, 
1 < m < d. The singular component of the distributions of X d (t) and Y d (t), t > 0, are projected 
on the subspaces R m , 1 < m < d, and enter into the absolutely continuous part of P{X.f n (t) £ 
dx m } and P{Yi(t) G dyj. 

Theorem 6. For t/ie random flights X m (t),t > 0, and Y m (t),t > 0, we have the following 
unconditional distributions 



P{Xl(t)edx m }_(c 2 t 2 -||x m || 2 )^ 



fej d — 1 d-r 



Mc 2 * 2 -||xJ| 2 )^- 
(ct)^- 1 ' 



n^Li dx 3 

with d > 2, ||x || < ct and 



irTf(ct) d - 2 



E d-i d (At) 



(3.13) 



P{Y.j(t) e dyj = {<?# H^JI 2 )^- 1 Vg V (<*)^ 



(3.14) 



with d > 3, llv II < ct. 

Proof. We observe that the projection of the uniform distribution on the surface of T-L^ t onto M d_1 

is obtained by means of the relationship dH^.sm.6 — YijZi d%j, where sin 6 = — c ^ dl — ■ 

Then, we obtain that 



fi d jxd-vt;0) 



r(f) 



ULl dx j 



{cty-^i ^c^-ita-il 



and by performing the following integrations 



V c2 * 2 -n?ui 2 f V c2t2 -\\^-2\\ 2 J r(^)( c 2 t 2 -iix ii 2 )^ 21 - 1 

dx m+1 --- fi (x rf 1 ,t;0)dx d ^ 1 ^-^ [ — ^^ , 

-VSHijp V« a -n»,- 3 ii' ^^ J ** (ct)^- 2 r(^) 



with ||x m || < ct, which corresponds to (2.261 and (2.27) for n = 0. Then, for the process X TO (t) 
we obtain that 

l\-j=i ax i n=0 

_ 1 1 y (At)" (c 2 t 2 -H^H 2 ) 1 ^^- 1 )-^ 

Trf B*-i a (At) ^, r(s±i(d - 1) + i^ 21 ) ( rt )(n+l)(d-l)-l 



1 (c 2 t 2 



r 



Bj d — 1 d - ?-, 



At(c 2 t 2 -||x m H 2 )^ 
(ct)< d - 1 ) 



(ct) 



d-2 



E d-l d (At) 



with ||x„J| < ct. 



Analogous considerations on the random flight Y m (t),i > 0, yield result (3.14| 



D 
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m = 1 



d= 2 



i £i i 

(c 2 t 2 -a; 2 )~2 2'2 



\(c z t z -ii ) 2 



El x (Ai) 



d = 3 



i E, 



( c 2 t 2 -||x 2 l| 2 r 2 ^H ^ 

-net E 1 3 (At) 



XtCeftf-llgallfj 






Table 4: Unconditional densities of X ? 2 „(t), t > 0, and X^((),i > 0. 
d=3 d=A 



to = 3 



to = 2 



TO = 1 



(c 2 * 2 -||y 2 ll 2 H ^ 



Xt(c^t^-||y ||2)2 



-El 



El 3 (At) 

2 ' 2 



A(c2 t 2_ H 2 )2 



^ B i,|(At) 



^2 + 2 II, , i|2\-i E 



At(c 2 t 2 -||y J| 2 )-s fi i,U (^ 



^At(<= 2 t 2 -||y 3 ll 2 )^ 



W5(ct) 2 



At ec2 t 



e A *-l 



(e-***-||sr,lr) 



7r(ct) 2 



e At -l 



w/ 2 + 2 2si E 3 

At(c t -j/ 1 ) 2 1, 2 



(ct)2 



wT (c t)2 - At -i 



Table 5: Unconditional densities of ~Y_ m (t),t > 0, and Y m (t),t > 0. 



(3.14) 



In Table 4][5 we sum up some particular important cases of the distributions (3.13) and 



Remark 3.3. The projection onto the one- dimensional space K 1 of (3.12) and of the singular 
component of a planar random flight with Poissonian times, becomes 



\e- xt ^ / A 



—2 / \ \~^ 



2c ^ V 2c v ' 

fc=0 



fc-i 



ra(*±l) 



c 



-At 



A 



e 



Try^F^? 2c £j 

e- xt Xer xt f /A 



e n V2cV V r^i + i 



V^T 2c ^o^v^-^J+^U^-^ 



(3.15) 



where \x\\ < ct, and Iq(x) = X^feLo fin 2 ' x ^ ^> * s ^ e "modified Bessel function, while Lo(x) 
XfcLo /p/ t .4_3-|-,2 i i£K, is the modified Struve function. It is particularly interesting to compi 
(3.15) with the probability distribution obtained in TableUlfor d — 2 and m = 1, namely 



P{Xl{t) € d Xl } 
dx\ 



Ex i 

2 '2 



(^^cW^xt) 



Wc 2 * 2 - x'f E hi( Xt ) V^cEi A (Xt) ~ V r 



E \-M^-x\ 



fc-1 



r(*±i) 



(3.16) 



2i : 



with \xi\ < ct. 

Remark 3.4. We observe that the k-th term, for k > 2, of (3.161 can be extracted from the 



uniform distribution inside the hypersphere H. ct — {x±, ...,Xk '■ | |x fc 1 1 < ct} as follows 



r(i) /V^* 2 -* 2 /•V c2 * 2 -H^-2ll 2 /V<= 2 * 2 -ll2£*-ill 2 
9k{x\) = — fc 2 / ^ cfe 2 --- / dajfc-i / dx k 

2ll 2 {ct) k J-y/cHl-xj J-y/cH*-\\x k _ 2 \\* J-V c2 * 2 -H2£ fc -lll 2 



r(| + i) 
Ar(^)(ct) fe 



2 £ 2 - x 2 



fc-1 



TTie above distribution for k — 2 yields the well-known Wigner law of which it represents an 
extension. We observe that 



/vFr(m +| + 1) 



ano 



a! /or fc = 2 yields J_ .x\ m g2{xi)dxi — ( m ) ,. ,2m+i 2(ct) 2 '"" 1 " 1 involving the Catalan num- 
bers and the distribution of the first return in the origin of the coin tossing process. 

By summing up the distribution gk{x\) with weighting terms represented by the fractional 
Poisson process ]\[<2{t),t > 0, with probability distribution 



we obtain the probability law (3.16). By integrating the uniform law in the hypersphere H^ t 
with respect to the variables x m+ i, ...,Xk, we obtain an m- dimensional extension of gk{x\) in the 
following manner 

r(fe) /Vc 2 * 2 -iix m ii 2 /V c2f2 -N**-2ii 2 /-V c2t2 -ii^-iii 2 

g fe (xi,...,a; m )= — F y — - / dx m+ i--- \ dxk-i / dx k 

2n 2 (ct) k J-y/w-w^w* •/- x / c2 * 2 -H^-2ll 2 J-V c2 * 2 -H^-ill 2 



r(f + i) 



k—m 



, V / c 2 t 2 -||x m || 2 ) , (3.17) 

^fr(^ + i)(ci) fc v v "-'"" ' 



with k > 77i. for m = fc — 1, (3.17) represents the m-dimensional extension of Wigner law. 



4 On three-dimensional random flights governed by a Pois- 
son process 

The space M 3 is an environment particularly important for the representation of the real mo- 
tions. Therefore, we will focus here our attention on the random flights developing in the three- 
dimensional Euclidean space. 

For the analysis developed in this Section, it is useful to observe that the Dirichlet distribution 



(1.4) is related to the Poisson process and thus permits us to give an interesting interpretation 
of the random flight. If Ti, ..., T n are the random instants at which the events of a homogeneous 
Poisson process occur it is well-known that 

P{7\ e dtu ...,T n G dt n \N{t) =n} = -dh ■ ■ ■ dt n l {0 < tl < t2 <...< 4 „< t} (4.1) 
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By integrating (4.1) as follows 



P{T d -i e dt d - 1 ,T 2 ( d - 1) e dt 2 ( d -i),-,T n ( d _ 1 ) € rf£»(rf-i)|JVft) = (n + l)(d- 1) - 1} 

dtd-\dt2(d-i) ■ ■ ■ dt n (d-i) 
_ ((n + l)(d -!)-!)! f 

~ t („+l)(d-l)-l / dtl ' ' ' "*«(d-l)+l " " " «*(n+l)(d-l)-l 



{0<ti<-<i,j_i<-<i n ( (1 _i)<-<t(„ + i) ((J _ !)_!<*} 

/ dti--- dt d ^ 2 / 



^ + D{d-l)-l)l ' -. > ' dt d -dt 



t (n+l)(d-l)-l _/ atl---«^-2 y «t«J---«l2(d-l)-l 

{0<t 1 <--<t d -2<td-l} {td— l<td<---<t2(d— 1) — l<*2(d-l)} 



/ 



rftn(d-l)+l ' ' ' ^(n+l)(d-l)-l 

{*n(d-l)<'n(d-l) + l<"'<*(n+l)(d-l)-l<'} 

{{n + l){d -!)-!)! tdli (t 2 ( d -i) - t d -i) d - 2 {t-t n{d _ 1) ) d - 2 
t ( n +i){d-i)-i (d-2)\ (d-2)! "' (rf-2)! 

r((n+l)(rf-l)) "+ 1 d _ 2 

(r(d- l))«+it(™+i)(rf-i)-i 11 i 



with Tj — ij(d-i) — i(j_i)(d_i), j = l,...,n + 1, we obtain the law (1.4). For some details on 
these calculations see Lachal et al. (2006). The derivation of the Dirichlet distribution given 
here permits us to describe the random flight as a motion in M. d where a Poisson process governs 
the changes of orientation every d — 1 events. In other words, every d — 1 events of the Poisson 
process the moving particle changes direction ignoring all the previous d — 2 events. For spaces 



of even dimension an analogous derivation of the second form (1.5) of the Dirichlet law can be 
envisaged. 

We introduce a random motion in M 3 slightly different from that introduced in the previous 
Section. We suppose that the changes of direction are governed by an homogenous Poisson 
process. In particular, we assume that the particle changes direction (uniformly distributed on 
the surface of the sphere) only at even- valued Poisson events. Therefore, if the number of Poisson 
events is N(t) = In + 1, n > 1, we have that the position of the particle at time t is represented 
by the following vector 

n+l n+l 



^3(*) = c^(t 2 k - t 2 k-2) sin6> 2 fc~2 sin0 2 fc-2 = c y~] r fc sin 9 2 k-2 sin 4>2k-2 
fc=i fc=i 

n+l n+l 

U 2 {t) = c^2(t 2k - t 2k - 2 ) sm6 2 k-2 cos 2 fc-2 =c^r t sin^-2COsfe_2 
fe=i fe=i 

n+l n+l 

Ul(t) = Cyifok - hk-2) COS 9 2k -2 = C^Tk COS 6 2k -2 



fe=l k=l 

where £& is the instant at which the fc-th Poisson event happens with t = and ^2fe+2 = £■ In 
view of the above considerations, for d = 3, one has that 

P{T 2 edt 2 ,T 4 edt 4 ,...,T 2n edt 2n \N{t) = 2n + l} _(2n + l)\ , £} 

dt2dU ... dhn - ^TT" 11 (*W - %-2) 

3=1 
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and then, the random flight U 3 (t) = (Ui(t),U2(i),Ua(t)),t > 0, has conditional characteristic 
function given by 

E je l< ^U3(t)>|AT(t) =2n+l\ 

/rj I i\| ft ft p27T p27T p7T p2tT 

= 2n+1 ' / t 2 dt 2 --- {t 2k - t 2k - 2 )(t - hk)dt 2 k d(f> --- 4>2n d9 --- 9 2n 

' Jo Jt 2k -2 Jo Jo Jo Jo 

ia 3 c^(i 2 fe - i 2 fe-2)sin02fc-2sin<?!>2fe-2 +«a2C^(t 2 fe - tik-2) sin0 2 fe-2 cos^ 2 fc-2 
fe=i fc=i 



n+l 



««lC^(/J2fc -t 2 k-2)cOs8 2 k-2 



fc=l 



sin #o sin 02 • • • sin t/2n 



(47T) 



n+l 



(by using (2.5) for d = 3) 



(2n + l)! /7T 



i 2 ™+! V2 



/ t 2 dt 2 (t 4 - t 2 )dt 4 ■ ■ ■ (t 2 k-t 2k - 2 )(t-t 2k )dt 2k Yl 

Jo Jt 2 Jt 2k -2 u--, 



n +}J h {c{t 2k -t 2k ^ 2 )\\a 3 \ 



{ ^^r\i;^sr— Ja 



2fc-2 

t-YZZl r k 



, =1 \/c(*2fc — *2fc-2)||^3l 



" n +l Ji(C7)fe||S3|r 

Tn(* - V r fc )dr„ I I 2 

fe =i fc =i V CT fcll« 3 ll 



(2n + 1)! /tt\ ^ j 2t+»+ 1 r(?i + |)i n+ 5 



T J„ + i(ci||a 3 | 



2«+sr(n+ |) 



^^J„+i(ci||a 3 | 



* 2 " +1 V2/ ^^r(2n + 2)( C ||a 3 ||)»+i " T ^ "-'"J (ci||« 3ll/ 

where in the last step we have used the same approach as that developed in the proof of the 



result (2.1 1. Hence, the characteristic function of U 3 (i),t > 0, (conditionally on the event 
N{t) — 2n + 1) coincides with ( |2.1[ ) for d = 3. Then, by inverting the Fourier transform 
E {e i< ^.u 3 (t)>| iV (^ _ 2n + i} ; we immediately obtain that 



P{U 3 (t) e du 3 |7Y(i) = 2n + 1} _ r(n + 1) 
rij=i^ 7r3r(n)(ci) : 



n+l v "— du ' ' — ' 



(4.2) 



with ||u 3 || < ci, and coincides with the distribution (2.101 for d = 3. For n = 1, the result (4.2) 



shows that we have an uniform distribution inside the sphere %%■ If n = 0, that is iV(t) = 1, 
the random flight U 3 (£),i > 0, (changing direction only at even- valued Poisson events) reaches 
the surface of %% with probability P{U 3 (i) £ dH 3 ct } = Xte xt . 

In this case, we can provide the following unconditional probability law 

P{V 3 (t)edu 3 ,l)Z 1 (N(t) = 2n + l)} = ~ p{N{t) = ^ + i] P{V 3 (t)edu 3 \N(t) = 2n+l} 

n5=i du J n=i ' n5=i dt *j 



-At °° 

E 



(Ai) 2n+1 T(n+|) 



.2+2 



tt! ^(2n + l)!I»(ci) : 



n+l v 



|2\n-l 



-At °° 



(At) 



2n+l 



/^r(2n + 2)2- 



-2n-l 



ttS ^(2n + l)!r(n + l)r(n)(ci) : 



n+l v ' ' — ■ i] ' ' 



2\n-l 



„-At 



7T V 2 C/ V^t 



U 



E 



{W< 



cH*-\\u 3 \ 
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r(n + l)r(n) 
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.,-At / A X 2 



/! -V^Hfe 



7T \2cJ V / C 2t2_||u 3 

Moreover, we obtain the projection of the absolutely continuous component of the distribution 
of U 3 (£),£ > 0, onto the plane as follows 

P{Ui{t) G dm, U 2 (t) € du 2 \N{t) = 2n + 1} _ /V c2 * 2 -I|u 2 || 2 P{U 3 (t) e du 3 \N(t) = 2n + 1} 
du 1 du 2 i_^/ c 2 t 2_||u 2 ||2 duidu 2 

n+ 4 T ( C 2 i 2 -||u 2 || 2 ) n -5 



7r(ci) 2 ™ +1 
for n > 1, while the projection of the uniform distribution on the surface &}i\ t onto the plane 

t*o *\ c\ ^ 

P{?7i(t) € dm, L^ 2 ft) g riM 2 |iV(£) = 1} _ 1 1 

duidu 2 2-irct y/cH 2 - ||u 2 || 2 ' 

Then 

P{Ui(t) € d Ul ,U 2 (t) € du 2 ,U~ (iV(*) = 2n + l)} 
duidu2 

= V P{7V(t) = 2n + 11 P{t/lW € ^ ^ € dualjVX*) = 2n + 1} 

n— 

e -At ~ (Ai) 2n + i n+ i N „ ||2 , , e"*< A~ 1 /AN 2 " 



7T 

n=0 



^ (^ixw^ '^ -""■"'»-' %^-iirf^S(^^ , ^ qHF|i 



cosh I — \/c?t 2 — ||u- H- 



2 CTv /c 2 t 2 -||u 2 || 2 ^""Vc V ^ "- 2 

We note that g = q{u u u 2 ,t) = P{U t (t) e du 1 ,U 2 (t) e dt^.lXLoW*) = 2n + 1 )} is a 
solution to the planar telegraph equation 

d 2 q „, <9o 9 f <9 2 <9 2 1 , , „, 

^ +2A J =c W^h (4 - 3) 

Furthermore, we have that 

P{Ux(t) e dui|AT(*) = 2n + 1} fV c2t2 - u2 i P{Ui{t) e du u U 2 {t) <= du 2 |JV(i) = 2n + 1} 



dtil 7_^ c 2 t 2_ u 2 C?Ul 

_ (2n + l)l 22 2 

~ (n!) 2 (2ci) 2 "+ l1 lJ 

and then 

P{E/i(t) ^ dm,U~ (iV(t) - 2n+ 1)} = P{T(t) 6 dm,U~ (JV(t) = 2n + 1)} 
du\ du\ 



(4.4) 



•^ l w J dwi ^ (n! 2 (2ci 2 "+ lV u 



= e -xt A 7 ("A / r2t2 _ u2 \ = i J {r(t)ed Ml ,U~ (JV(t) = 2n + i)} 

2c ° V2cV V dm 
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The result (4.4 1 shows that the projection of U 3 (t),t > 0, onto the real line is equivalent in 
distribution to the classical telegraph process T(t),t > 0. 

If the Poisson events recorded in [0, t] are 2n, we are not able to express in closed-form 
P{U 3 (t) G du 3 \N(t) = 2n}. Indeed, we have that 

P{U 3 (t) G du 3 \N(t) = 2n} 

n5=i du i 



yXx [ e - l< ^^3> E \e l< ^^ it)> \N(t) =2n\f\da j 



Hence, to evaluate explicitely the distribution -P{U 3 (£) G du3|iV(i) = 2n}, we need to calculate 
the integrals of the following form 



m+r 2m+2r+2 

yfa-^c"^'"^ "'"" ~ ^ ^ m\r\T(m + \ + l)r(r + \ + 1) 2 2 -+ 2 ™+i ,,„ 



J, (x) J, (a - x)dx = E E „, Mr L fr-1 u^ rr-T^ 92r+2m+1 / l^Ha-y)^ 



oo oo 



EE 



(-l) m+r a 2m+2r + 2 1 r(2m + 2)r(2r + l) 

=o ^o w!r!r(m+| + l)r(r+| + l) 2 2 ''+ 2 ™+ 1 r(2(m + r) + 3) 

o °° °° (" i\r+m 2(m+r)+2 



77 ^-^(2r + l)r(2(m + r) + 3) 

m=0 r— ' 

~ 7T ^ ^ (2(r - m) + l)r(2r + 3) ~ tt ^ T(2r + 3) ^ 2m + 1 

m— r—m v v / / \ / r _g \ / m _ q 

Unluckily, the above integral cannot be worked out explicitly and then the recursive approach 
used in the proof of Theorem fll cannot be applied. Nevertheless, if N(t) = 2, since Ji(x) — 

— sinx, we have that 

P{U,(t) G du 3 \N(t) = 2} = du ^du 3 1 r« Jr^ r~ sinCHIualDsinCcr^BinCcft-rOp) 

tt 2 (c£) 2 ||u 3 || J a t- Tl J p 

The integral with respect to p can be treated as in Orsingher and De Gregorio (2007). Indeed, 
since 

sin x sin y sin z — — [sin(z + x — y) + sin(z — x + y) — sin(z + x + y) — sin(z — x — y)\ 



and in light of the remarkable fact that 

An 

-dp 



sia(Ap) j_ J f, A>0, 



■f, ^<0, 
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the integral becomes 



1 J -^| sin(p(ct - HU3ID) + sin(p(c* - 2csi + H1I3II)) - sin(A>(ct + HU3ID) - sin(p(ct - 2e»i - HU3II))} 



2 2 7o 



7T 

= ^-lr ct-ii„3ii^f+ii„ 3 in (ti) 

Then 

P { U 3 (i) e d ^N(t) = 2} = %^^ log (|±|g) , llu.ll < ct, 

which corresponds to the law obtained in IR 3 , for n = 1, in Orsingher and De Gregorio (2007). 
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